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Abstract. In strongly correlated materials the electronic and optical properties are 
significantly affected by the coupling of fermionic quasiparticles to different degrees of freedom, 
such as lattice vibrations and bosonic excitations of electronic origin. Broadband ultrafast 
spectroscopy [I] [2] is emerging as the premier technique to unravel the subtle interplay 
between quasiparticles and electronic or phononic collective excitations, by their different 
characteristic timescales and spectral responses. By investigating the femtosecond dynamics 
of the optical properties of Bi2Sr2Cao.92Yo.osCu208+5 (Y-Bi2212) crystals over the 0.5-2 eV 
energy range, we disentangle the electronic and phononic contributions to the generalized 
electron-boson Eliashberg function 3 4, showing that the spectral distribution of the electronic 
excitations, such as spin fluctuations and current loops, and the strength of their interaction 
with quasiparticles can account for the high critical temperature of the superconducting phase 
transition [5]. Finally, we discuss how the use of this technique can be extended to the 
underdoped region of the phase diagram of cuprates, in which a pseudogap in the quasiparticle 
density of states opens. 

The microscopic modeling of the interaction of ultrashort light pulses with unconventional 
superconductors will be one of the key challenges of the next-years materials science, eventually 
leading to the full understanding of the role of the electronic correlations in controlling the 
dynamics on the femtosecond timescale. 



1. Introduction 

The generalized bosonic function 

In conventional metals the scattering between quasiparticles (fermions) and lattice vibrations 
(bosons) is the microscopic mechanism that determines the transport and the optical properties. 
At low temperatures the electron-phonon interaction, denoted by a 2 F(Q), provides the "glue" 
for the formation of the Cooper pairs, leading to the instability of the Fermi-liquid ground 
state upon the formation of the superconducting condensate. The critical temperature of 
the superconducting phase transition, that is of the order of a few degrees, is determined by 
the electron-phonon coupling constant, i.e. A e _ p ^=2 f a 2 F(Q)/Q d£l, through the McMillan's 
formula [6]. 

In the cuprate superconductors, the electronic correlations are responsible for the emergence 
of complex orders, influencing the phase diagram even far from the antiferromagnetic insulating 
phase at zero doping. These novel degrees of freedom, like paramagnon excitations [7], charge 
fluctuations [8] or loop currents [9] provide additional scattering channels to quasiparticles 
(QPs), strongly affecting their lifetime and dispersion over an energy range of the same 
order of the spectral distribution of the excitations. An "effective" description of the generic 
interaction between quasiparticles and excitations of both electronic and phononic nature, whose 
distribution at a given temperature follows the Bose-Einstein statistics, can be obtained by 
replacing the a 2 F(Q) with a more general electron-boson coupling function i.e., the Bosonic 
Function 11(0), defined as: 

u(n)^a 2 F(n) + i 2 x (n) (i) 

where l 2 x(0) accounts for the coupling with all the bosonic excitations of electronic origin. 

The signatures of the QP-boson interaction usually manifest in most of the experiments that 
probe the electronic properties at equilibrium. The generalized bosonic function II(O) can be 
extracted [I] by analyzing the kinks in angle-resolved photoemission data [10] , the dip features in 
tunneling spectra [HJ Q2] , the frequency-dependent scattering rate in optical spectroscopy [131 E] 
and Raman spectra [13]. 11(0) shows some ubiquitous features, independent of the technique 
used, like a strong peak at 50-80 meV and a broad continuum that extends up to 300-400 meV. 
Although a main effort has been recently made [1] , the determination of the relative weight of the 
electronic and phononic contributions to 11(0) remains elusive, since the electronic and phononic 
excitations coexist on the 0-90 meV energy scale. Solving this problem would constitute a step 
forward to unravel the puzzle of superconductivity in cuprates, while addressing the major 
question whether or not high-temperature superconductivity can be described in terms of a 
generalized Eliashberg formalism, in which the attractive interaction is 'retarded', i.e., mediated 
by virtual bosonic excitations of novel origin |15] . 

Non-equilibrium optical spectroscopy is emerging as a new tool to tackle these fundamental 
questions [U [2] . The founding concepts of this technique rely on the use of an ultrashort light 
pulse (pump) to prepare the system in a non-equilibrium state, i.e., with the distribution of 
fermionic QPs decoupled from the distribution of bosonic excitations. The induced transient 
changes in the optical properties can be modeled by a two-steps process: i) the sudden photo- 
injection of fermionic QPs results in an effective increase of the plasma frequency of the Drude 
peak in the optical conductivity, without any change of the QPs scattering rate; ii) the subsequent 
heating of the bosonic excitations induces an increase of the scattering rate and a broadening of 
the Drude peak, without any change of the plasma frequency. 

Since changes in the plasma frequency and in the scattering rate affect the optical conductivity 
in qualitatively different ways, optical spectroscopy with femtosecond time-resolution can be used 
to probe the temporal evolution of the distribution of fermionic QPs and bosonic excitations, as 
a function of the delay from the pump pulse. Furthermore, the dynamics of bosonic excitations 
of electronic and phononic nature can be disentangled on the basis of their different timescales 
while exchanging energy with QPs. 
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Figure 1. Reflectivity of the Bi2Sr2Cao.92Yo.o8Cu208+<5 crystal at optimal doping (T c =96 
K), as measured by conventional spectroscopic ellipsometry (red dots). The reflectivity is well 
reproduced by an extended Drude model and a sum of Lorentz oscillators (black line) . The total 
bosonic glue, n(O), extracted from the data is shown in the left inset. The dashed black line 
represents the reflectivity calculated for the same system at T=500 K. 

In this work, we will briefly introduce the basics of equilibrium and non-equilibrium optical 
spectroscopy, we will review the main results on the determination of the relative electronic and 
phononic contributions to 11(0,) and, finally, we will show how this technique can be extended 
to investigate the opening of a pseudogap in the quasiparticle density of states and possible 
anomalies in the temperature dependance of the QP self-energy. 

2. Optical spectroscopy at equilibrium 

Optical spectroscopy at equilibrium is a fundamental tool to investigate the electronic properties 
of strongly correlated materials |16j and unconventional superconductors |17j . Since the dielectric 
function contains direct information about the scattering rate of QPs and their dynamical 
effective mass, optical spectroscopies are intrinsically sensitive to the coupling of QPs to bosonic 
excitations. 

2.1. The optical properties of a prototypical cuprate 

In Figure 1 we report the aft-plane reflectivity of optimally-doped Bi2Sr2Cao.92Yo.o8Cu208+<5 
(Y-Bi2212) high-quality crystals |18j (T c =96 K), measured by conventional spectroscopic 
ellipsometry |19| at 300 K. The normal-incidence reflectivity R(uj, T) is related to the dielectric 
function e(u,T) by the relation: 



R(lo,T) 



(2) 



The dielectric function is obtained by applying the Kramers-Kronig relations to the reflectivity 
for 6 meV<^o;<0.74 eV and directly from ellipsometry for 0.2 eV</iu;<4.5 eV. This combination 
allows a very accurate determination of e(oo) in the entire combined frequency range. 



The reflectivity of Y-Bi2212 shows some general features common to most of cuprates: 

i) below 1.25 eV the optical properties are dominated by a broad peak related to the optical 
response of the low-energy excitations in the conduction band. Interestingly, the dielectric 
function in this energy range cannot be simply reproduced by a Drude peak, in which a constant 
scattering rate r and effective mass m* of the QPs is assumed. The strongly frequency-dependent 
scattering rate and effective mass, resulting from the interaction with bosonic excitations, will 
be accounted for by the more general extended Drude model |17j . presented in the next section; 

ii) above 1.25 eV, the high-energy interband transitions dominate. In this energy range, the 
equilibrium dielectric function can be modeled as a sum of Lorentz oscillators at ~1.5, 2, 2.7 
and 3.9 eV [2]. The attribution of these interband transitions in cuprates is a subject of intense 
debate. The ubiquitous charge-transfer (CT) gap edge (hole from the upper Hubbard band with 
d x 2_ y 2 symmetry to the 0-2p x>y orbitals) in the undoped compounds is about 2 eV [IT]. Upon 
doping, a structure reminiscent of the CT gap moves to higher energies, while the gap is filled 
with new transitions. This trend has been recently reproduced by Dynamical Mean Field Theory 
(DMFT) calculations of the electron spectral function and of the afr-plane optical conductivity 
for the hole-doped three-band Hubbard model [20]. The structures appearing in the dielectric 
function at 1-2 eV, that is, below the remnant of the CT gap at 2.5-3 eV, are possibly related 
to transitions between many-body Cu-0 states at binding energies as high as 2 eV (for example 
singlet states) and states at the Fermi energy; 

hi) the dressed plasma frequency, uJ p , defined through the relation Re{e(cJ p , T)}=0, is ~1 eV. In 
this energy range, the signatures of the pump-induced change of the QP distribution and of the 
heating of bosons can be more effectively disentangled, as will be discussed in Section 3. 



2.2. The extended Drude model 

In the extended Drude model (EDM) the physical processes responsible for the renormalization 
of the lifetime and effective mass of the QPs are accounted for in a phenomenological way, 
by replacing the frequency-independent scattering time r with a complex temperature- and 
frequency-dependent scattering time r(u,T): 

t- 1 => f ~V) = t~ x (oj) - *wA(w) = -iM(oo, T) (3) 

where 1 + A(w) = ^r(w) is the mass renormalization of the QPs due to many-body interactions 

and M(uj,T) is the memory function. 

In the EDM, the dielectric function €d(uj,T) is given by: 

e D (u,T) = 1 - "p 2 , . . = 1 (4) 

while the Drude optical conductivity crjj(oj,T)=[l — €d(uj, T)]iuj/4ir reads: 

cr ( u T)= i Up2 = — ^ (5) 



The renormalized scattering rate and effective mass can be directly extracted from the measured 
Drude optical conductivity, through the relations: 

1 + A(„,T) _ -^I, m (_±_) (7) 




Figure 2. Frequency-dependent 
scattering rate (left axis) at 300 K 
(red line) and 100 K (blue line) 
and effective mass (right axis) at 
300 K (red dashed line) and 100 K 
(blue dashed line) of the fermionic 
quasiparticles. 
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Although this phenomenological version of the EDM does not provide any clue about the 
microscopic mechanisms responsible for the renormalization of the energy dispersion and lifetime 
of the QPs, it is very useful for directly extracting t(uj, T) and m*(uj, T)/m from the optical data. 
In Figure 2, we report t(uj,T) (solid lines, left axis) and m* (u,T)/m (dashed lines, right axis) 
obtained by applying Eqs. 6-7 to the low-energy part of <7d(w, T), measured on optimally-doped 
Y-Bi2212 samples. At T=300 K (red lines), the value of the QP scattering time is ~5 fs at very 
low frequencies and decreases towards an asymptotic value of ~2 fs above 0.5 eV. These values 
of t(u, T) imply that QPs scatters very quickly and exchange energy with bosons on the very 
femtosecond timescale. As the temperature decreases, the opening of the pseudogap prevents 
QPs from being scattered and an increase of the scattering time below 100 meV is measured, 
while the value of t(ui,T) above 100 meV is almost unaltered. A similar behavior is measured 
for the effective mass, that decreases towards the m*/m=l asymptotic value above 0.5 eV. 



2. 3. The extraction of the bosonic function from the optical conductivity 

From the microscopic point of view, the Extended Drude formalism can be derived from the 
Holstein theory for normal metals |21j . Considering the Kubo formula and using complex 
diagrammatic techniques to evaluate the electron and boson thermal Green's functions and 
omitting vertex corrections (Migdal approximation), the Memory function M(oj,T) results: 



M(uj,T) 



UJ 



m,T)-m+Lj,T) 

+ £*(£,T)-£(£ + W ,T) + 2 7imp 



d* 



-i 



— UJ 



(8) 



where / is the Fermi-Dirac distribution, E(cj,T) and S*(w,T) the electron and hole k-space 
averaged self-energies and 7i mp an intrinsic decay rate that accounts for the scattering by 
impurities. We pinpoint that, although the memory function M(uj,T) has the same analytical 
properties of the single-particle self-energy S(w,T), it has a conceptually different meaning, 
since the optical transition at frequency uj involves a particle-hole excitation of the many-body 
system and provides information about the joint particle-hole density of states. 

The electron self-energy T,(uj,T) can be calculated as a convolution integral between the 
bosonic function n(O) and a kernel function L(uj,£1,T): 



The kernel function 
L(u,n,T) 



s(w, t)= I u(n)L(u),n,T)dn 

'o 



ra(n',T) + /(n,T) l + n(Q',T) - f(Q,T) 



Q - uj + VL' + i5 



+ 



Q — uj — Q,' — id 



(9) 



(10) 



accounts for the distribution of the bosonic excitations through the Bose-Einstein distribution 
n(0,T), and can be calculated analytically: 



L(lu, 0; T e 



T 6 ) = -2vri n(0,T b ) + J 



O — oo 



2irT e 



) 




2 2vrT e J V ' 



where \P are digamma functions and the dependence of the different terms on the temperatures 
of the electronic QPs (T e ) and bosonic excitations (T&) has been made explicit. In this formalism, 
the frequency-dependent scattering rate is a consequence of the microscopic interaction of QPs 
with a distribution of bosons at temperature T&. 

Close to optimal doping, the vertex corrections beyond Eliashberg theory can be reliably 
neglected and the EDM can be safely used to extract II(O) from the optical conductivity, 
measured at the equilibrium temperature T. Although sophisticated maximum entropy 
techniques p3] have been developed to unveil the rich details of the bosonic function, the main 
features can be evidenced by simply fitting the model in (5) to the experimental reflectivity 
R(ui,T), reported in Figure 1, and assuming a simple histogram form for 11(0). The resulting 
bosonic function is reported in the left inset of Figure 1 and is characterized by: 

i) a low-energy part (up to 40 meV) linearly increasing with the frequency. This part is 
compatible with either the coupling of QPs to acoustic [22] and Raman-active optical |23| 
phonons or the linear susceptibility expected for a Fermi liquid |24| : 

ii) a narrow, intense peak centered at ~60 meV, attributed to the anisotropic coupling to either 
out-of-plane buckling and in-plane breathing Cu-0 optical modes [25J or bosonic excitations of 
electronic origin such as spin fluctuations [26j; 

hi) a broad continuum extending up to 350 meV, well above the characteristic phonon cutoff 
frequency (~90 meV) and usually attributed to the coupling with spin fluctuations [271 [28], \29 \ IT] 
or loop currents [30] , 

Using a sum of the EDM and four Lorentz oscillators at ~1.5, 2, 2.7 and 3.9 eV, the reflectivity 
of optimally-doped Y-Bi2212 can be satisfactorily reproduced up to 2.5 eV photon energy (see 
the black solid line in Figure 1). In the EDM the main role of the temperature is to change the 
density of the bosons and, as a consequence, the scattering rate of QPs. In Figure 1 we report 
R(oo,T) calculated for T=500 K (dashed line). The T-related increase of the scattering rate 
induces a broadening of the Drude peak, resulting in a decrease of the reflectivity below oo p =l 
eV and an increase of the reflectivity above oo p . 

Although very useful to determine the microscopic origin of the frequency-dependent 
electronic properties and capable of reproducing most of the features of the temperature- 
dependance of the dielectric function, the EDM does not allow to evaluate the relative 
contributions of the phononic and electronic contributions to the total II(u;), since they spectrally 
coexist below 90 meV. 

3. Non-equilibrium optical spectroscopy 

Non-equilibrium optical spectroscopy permits to monitor the change of the dielectric function 
and, in particular, of the QPs scattering rate, during the relaxation process after the interaction 
with an ultrashort optical light pulse. Adding the temporal (t) dimension to the frequency (oo) 
dimension will be the key to disentangle the phononic and electronic contributions to H(cj). A 
sketch of the technique is reported in Figure 3. A 100 fs infrared optical pulse (1.5 eV photon 
energy) is used to excite the sample. The dynamics of the dielectric function is probed by 
combining two ultrafast optical techniques: i) supercontinuum light generation in a photonic 
fiber and simultaneous detection of the 1-2 eV spectral region through a linear array, after 
dispersion in a prism; ii) Optical Parametric Amplification providing tunable and short (<100 
fs) output pulses in the 0.5-1.1 eV range. The details of the experiment can be found in Refs. 




Figure 3. Sketch of the experimental 
setup for the 1.5 eV-pump and 
colored-probe optical measurements. 



3. 1 . The extended Drude model out of equilibrium 

During the interaction with the pump pulse, two physical processes are expected to take place 
on the sub-picosecond timescale: 

-High-energy fermionic excitations are impulsively injected by the pump pulse. Since the lifetime 
of these excitations is very short, they rapidly relax toward the Fermi energy, producing a high- 
density of QPs at lower excitation energies. As far as no pseudogap is present, which would 
provide constraints for the QPs scattering processes, we assume that the excess of QPs can be 
described through an effective increase 5T e of the equilibrium electronic temperature T e . 
-The temperature T b of the distribution of the bosonic excitations shows a 5T b increase because 
of the QP-boson coupling. 

By noting that the Kernel function L(oj, O,; T e , T b ) is the sum of separate terms depending on 
the electronic and bosonic temperatures, it is possible to use the EDM described in Section 
2.3 to calculate the changes in the optical properties expected at i=0, i.e., soon after the 
excitation, in the two limits: a) Non-thermal scenario, i.e., 5T e ^>5T b , in which the QPs are 
almost completely decoupled from the bosons, b) Quasi-thermal scenario, i.e., 5T e ~5T b , in 
which the QPs and bosons are in thermal equilibrium at an effective temperature larger than 
the equilibrium temperature. 

To facilitate the comparison with the experimental outcomes, we calculated the relative 
variation of the reflectivity at t=0, defined as the relative difference of the reflectivities out 
of equilibrium (R oe q) and at equilibrium (R eq ): 

SR. R oeq (uj, T) - R eq (uj, T) = R(cj, T e + 5T e , T b + 5T b ) - R{oo, T e , T b ) 

R [U ' ) Re q (u,T) R(u:,T e ,T b ) [ > 

The results are reported in Figure 4. In the case dTe^STf, (yellow solid line), the variation of 
the reflectivity shows a positive, intense feature just below uJ p =l eV, while it is negligible at 
higher energies. This behavior can be rationalized in terms of a small increase of the plasma 
frequency without any change in the scattering rate t(uS). For the case 5T e ~5T b (black solid 
line), the reflectivity variation is dramatically different, since it is dominated by the increase of 
the density of the bosons and of the QPs scattering rate, as evidenced by the broadening of the 
Drude peak. In the quasi-thermal scenario, 5R/ R(uj,T) turns from positive to negative at uJ p . 

These results open the way for probing the dynamics of the dielectric function over a wide 
energy range, above and below uJ p , allowing to discriminate between the non-thermal (5T e ^$>5T b ) 
and quasi-thermal (5T e ~5T b ) scenarios and directly measure the time necessary for the heating 
of the bosonic excitations. 

3.2. The bosonic function and the relaxation dynamics 

The key point to analyze the time-resolved data is that the same bosonic function Ii{Q) 
obtained from the extended Drude model (see Eq. 9) controls the temporal dynamics of 
the energy exchange between QPs and bosons, as pointed out by P. Allen in Ref. |31j . 



This model, introduced for simple metals and referenced as the "effective-temperature model", 
can be extended to capture the more complex physics of cuprates, including the strong 
coupling with some Cu-0 optical phonons. Formally, the total bosonic function can be written 
as U(£l)=I 2 x(Q)+ct 2 F(Q)scp+® 2 F(Q)iat, where I 2 x(ty refers to the bosonic excitations of 
electronic origin at the effective temperature T^ e , a 2 F(Q)scP to the small fraction of strongly- 
coupled phonons (SCPs), e.g., buckling and breathing Cu-0 optical modes, at Tscp[32] and 
a 2 F(yL)i at to all other lattice vibrations at T[ at , including acoustic and Raman-active optical 
phonons. 

A set of four coupled differential equations can be used to represent the following physical 
processes: a short laser pulse, with power density (absorbed) p, impulsively increases the effective 
electronic temperature of the QPs with a specific heat C e =7 e T e ('y e =ir 2 N c N(eF)k 2 /3, N c being 
the number of cells in the sample and N(Ep) the density of states of both spins per unit cell). 
T e will then relax through the energy exchange with all the coupled degrees of freedom that 
linearly contribute to the total U(Q)=I 2 x{Q)+ot 2 F(Q)scp- l rC i2 F(£l)iat- The rate of the energy 
exchange among the different populations is given by [31j : 

8T e _ G(I 2 X ,T be ,T e ) | G(a 2 Fscp,T SC P,T e ) | G(a 2 F lat , T lat , T e ) | p 



dt le T e le T e le T e le T ( 

dT be G(I 2 X ,T be ,T e ) 
8t C be 

dTscp G(a 2 Fscp,T S cp,T e 



dt Cscp 
dT lat G(a 2 F lat ,T lat ,T e 



dt Ci at 

where 



(14) 
(15) 
(16) 



G(U h Ti,T e ) = / dnUi(n)n 2 [n(n,Ti) - n(Q,T e )] (17) 

irfUtf Jo 

with Ui=I 2 x, a 2 Fscp, a 2 Fi at and n(Q, Tj)=(e n / fcsTi — 1) _1 the Bose-Einstein distribution at 
the temperatures Ti=T be , T$cp, Tlat- The specific heat {Cscp) of SCPs is proportional to their 
density of states and is taken as a fraction / of the total specific heat, i.e., CscP=fC[ a t- 

Combining the four temperature model (4TM, Eqs. 13-17) with the extended Drude model 
(Eqs. 5,8,9,11), provides a full picture of the dynamics of 8R/ R(co,T be ,Tscp,Ti at ), in which 
T be (t), Tscp(t) and Ti at (t) are functions of the time variable t and their dynamics is controlled 
by the three different ratios G(IIi, T|, T e )/Cj, where Ci=Cb e , Cscp and C\ a f 



4. Disentangling the electronic and phononic contributions to the bosonic function 

Figure 4 reports the SR/R(uj,t), as measured through non-equilibrium spectroscopy on 
optimally-doped Y-Bi2212 samples at T=300 K, versus time (i, x-axis) and frequency (uj, y- 
axis). Some worthwhile features emerge from the data. 

Already on the very short timescale (<100 fs), 5R/R(u,t) is negative below 1 eV and positive 
above 1 eV. The red dots are the values of 5R/R(u,t = 0), i.e., the maximum reflectivity 
variation versus the probe wavelength. The data are well reproduced by the reflectivity variation 
calculated in quasi-thermal conditions (5T e ~5T b ) through the EDM. It is worth noting that 
no signature of a change in the plasma frequency, characteristic of the non-thermal scenario 
(5T e ~^>5T b ), is detected within the time resolution of the experiment. Hence, it is possible to 
conclude that, already on the 100 fs timescale, the QPs are thermalized with some bosonic 
excitations participating to n(J7). The fast timescale of this thermalization implies a very large 
coupling and a relatively small specific heat. This finding strongly suggests that this process 




Figure 4. Relative reflectivity variation, i.e. 5R/ R(u,t)=(R exc (u),t)-R eq (uj))/ R eq (uj), as a 
function of the probe photon energy and delay between the pump and probe pulses. The data 
have been taken on an optimally-doped Y-Bi2212 sample at room temperature. For further 
details, see Ref. [5]. The maximum 5R/ R(u,t) at t=0 is calculated in the quasi-thermal (black 
line; T e ~T;,) and non- thermal (yellow line; T e 3>Tfe) scenarios, using the parameters obtained 
from the fit to the equilibrium measurements. The red dots are the 5R/R(oj, t=0) measured by 
time-resolved optical spectroscopy. 



involves bosonic excitations of electronic origin and is consistent with the scattering time t(cl>)~1- 
5 fs estimated from the equilibrium optical conductivity measurements reported in Figure 2. 
The temporal dynamics at fixed wavelength, i.e., 5R/R(t), is characterized by two decay times 
(~200 fs and ~1 ps), universally measured in time-resolved experiments [33, 34J. The faster 
dynamics is related to the coupling of QPs to SCPs, controlled by the a 2 Fscp part of the 
bosonic function and by Cscp, while the slower dynamics is attributed to the coupling of QPs 
to all other lattice modes and it is controlled by a 2 Fi at and Ci at . 

A quantitative analysis of the data can be carried out by using the combination 
of the four temperature model (4TM, Eqs. 13-18) and of the extended Drude model 
(Eqs. 5,8,9,11), previously discussed. The different contributions to n(S7) are extracted 
by fitting the calculated 5R/R(uj,t) to the time- and frequency-resolved reflectivity data, 
reported in Figure 4. Considering that the energy distribution of phonons is limited to 
<90 meV, we assume that, for f2>90 meV, II(0)~/ 2 x(^)- Within this assumption, the 
functional dependence of 5R/ R(u,t) on I 2 x(Q), a 2 F(Q,)scP and a 2 F(Q,)[ at is simplified as a 
parametric dependence on three coefficients pi, where I 2 x(^)=Pin(S7<90 meV)+Il(O>90 meV), 
a 2 -F(f2)scP=P2n(O<90 meV) and a 2 F(^l)i at =psIl(Q<90 meV). Considering the constraint 
that Tl(fl)=I 2 x{Q)+a 2 F(Q,)scP+a 2 F{tyiat (n(f2) being the total glue function extracted from 
equilibrium optical spectroscopy and reported in Figure 1) and fixing the values C e /T e =7 e =lCT 4 
Jcm~ 3 K -2 and Ci a t=2.27 Jcm _3 K _1 , the parameters phase-space of the model is significantly 
narrowed, allowing to unambiguously haul out the different contributions to 11(17) and to 



estimate Cfe e and Cscp- For further details, see Ref. [5]. 

The analysis of the 5R/ R(uj,t) demonstrates that the entire high-energy part and ~46% of 
the peak at ~60 meV instantaneously thermalize with the QPs at a temperature 2b e — T e . The 
spectral distribution and the estimated value of the specific heat of these excitations (C^ e <0.1C e ) 
support their electronic origin. The a 2 F(Q)scP is estimated to be ~34% of the peak at ~60 
meV, corresponding to ~17% of the total 11(f)), while a 2 F(Q)i at provides ~20% of the peak. 

The measured values of the QPs-boson couplings can be used to estimate the upper bound 
for T c , expected for each subset of the bosonic function. The maximal critical temperatures 
attainable are calculated assuming that ILj(Sl) entirely contributes to the d-wave pairing and 
fixing the pseudopotential fj,*=0. In the strong-coupling formalism, the critical temperature for 
d-wave pairing in a Fermi liquid with /i*=0, is approximately given by [6]: 

T c = 0.83O exp[-1.04(l + X l )/gX i ] (18) 

where lnfi=2/A; / °° II 6 (fi)lnft/JM), A;=2 / IT (ft) /0 dO is the electron-b oson coupling constant 
and ge[0,l] is a parameter that accounts for the d-wave nature of the superconducting gap. The 
upper bound g=l is used in the following. 
We summarize the main results: 

-The coupling to bosons of electronic origin is A e _?, e =2 J I 2 x(ft)/f2 d0=l.l±0.2. This value 
corresponds to a maximal T c =105-135 K. 

-The coupling to strongly-coupled phonons (most likely breathing and buckling Cu-0 optical 
modes) is A e _scp=2 J cx 2 F(£1)scp/Q dft=0.4±0.2. This value is in complete agreement with 
the values measured on similar materials via different techniques, such as time-resolved photoe- 
mission spectroscopy [32) . time-resolved electron diffraction [35] and single-color high-resolution 
time- resolved reflectivity [36J. The maximal critical temperature estimated considering only 
SCPs as the glue is T c =2-30 K. 

-The coupling to all other phonons is A e _; a j=2 J a 2 F(fl)i at /Q, d$7=0.2±0.2, corresponding to 
T c =0-12 K. 

The large error bars in the values of \ e -bei M-scp and M-iat include: i) the experimental 
uncertainty in the pump fluence and other experimental parameters; ii) the possibility of adding 
to the 4TM a term that accounts for the anharmonic coupling of SCP to the lattice; hi) the 
possible overestimation of Xscp an d Mat by a factor 8/5 [37], in the extreme case that the 
electron-electron scattering time (T e _ e ) is much larger than the electron-phonon scattering time 
and the energy exchange between the QPs and bosons begins before the establishment of a 
quasi-equilibrium fermionic population at T e + 5T e . 

Non-equilibrium optical spectroscopy allows disentangling the electronic and phononic 
contributions to the total bosonic function 11(17) of a cuprate superconductor. The strength 
of the coupling to bosons of electronic origin and its spectral distribution can account alone 
for the high-critical temperature (T c =96 K) of the system. This supports the description of 
high-temperature superconductivity in cuprates in terms of a generalized Eliashberg formalism, 
in which the attractive interaction is mediated by virtual bosonic excitations of electronic origin, 
such as spin fluctuations or current loops. 

5. Opening of a pseudogap and failure of the extended Drude model 

In the calculation of the self-energy in Section 2.3, a constant density of states at the Fermi 
level has been considered. Although this approximation is reliable at T=300 K in optimally 
and overdoped systems, it fails dramatically as the temperature and the doping decrease and a 
pseudogap opens in the electronic density of states. A further evolution of the EDM, accounting 
for a non-constant electronic density of states, has been developed by Sharapov and Carbotte 
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Figure 5. 5R/R(oj,t=0) calculated in the case of: impulsive heating of the fermonic 
quasiparticles and bosonic excitations (black line), impulsive filling of the pseudogap (blue 
line), impulsive quench of the electron-boson coupling (red line). The right insets qualitatively 
show the change of the density of states, N(u), and of the total glue, II(f2), assumed in the 
calculation of SR/R(uj, t=0). The equilibrium temperature Tn=100 K and the same n(S7) as the 
one determined at T=300 K have been used. 



[38], and has been used to analyze spectroscopic data at equilibrium [39 . In this model, the 
imaginary part of the electronic self energy is given by: 



(19) 



imS(oj, t) = -7T / n(n){N(uj + n, t) [n(n, t) + /(w + n, t)] + 

Jo 

+ n(u - n, t) [i + n (n, t) - f(u - n, T)]}dn 



while Re£(u;, T) can be calculated through the Kramers-Kronig relations. To mimic the opening 
of a gap, the normalized density of states N(oj,T) can be written in the following form |39j : 



N(Q,T) + [l-N(p,T)] 
l + |[l-iV(0,T)] 



for \u\ ^ A pg 
for\uj\ £ (A pg ,2A pg ) 
for 2 A pg 



(20) 



Where A pg is the energy gap width, while iV(0, T) is the gap filling. In the top-right inset of 
Figure 5 we report N(u,T) for the values A pg =40 meV and iV(0,T)=0.7. 

Replacing Eqs. 9 and 11 with Eqs. 19 and 20 it is possible to calculate the 5R/R(uj,t) 
induced by the pump pulse, even in the presence of a pseudogap. 
In particular, three different relevant processes can be analyzed: 

a) the impulsive heating of bosons at a temperature higher than the equilibrium temperature 
To of the system. A quasi-thermal scenario, with T e ~T{,, and Tn=100 K is assumed. The 
calculated 5R/R(oj,t = 0), reported in Figure 5 (black line), is very similar to the reflectivity 
variation obtained in the case of a constant density of states at T=300 K (see Figure 4) and is 
related to the impulsive increase of the scattering rate and broadening of the Drude peak; 

b) the impulsive filling of the pseudogap, as a consequence of the photoinjection of excess 
excitations. The 5R/R(u,t = 0) is calculated assuming that the value iV(0,T) is impulsively 
quenched from 0.7 to 0.9. The corresponding change in the density of states is shown in the 
top-right inset of Fig. 5. The result (see blue curve in Figure 5) is very different from the 



5R/R(uj, t = 0) expected for an impulsive heating of the bosons. Very similar results are obtained 
by assuming the closing of the gap (A p9 <A p9 o) instead of its filling; 

c) the impulsive decrease of the total II(J1). Although a 2 F(ft)scP and a 2 F{Q)i at are expected 
to be doping- and temperature- independent, I 2 xiS^) could increase as the doping and the 
temperature decrease because either new magnetic excitations coupled with QPs emerge in 
the pseudogap phase HQ] or electronic correlations induce a non-Fermi liquid like increase 
of the self-energy at the antinodes. 5R/R(u},t = 0) is calculated by assuming a small quench 
of the peak at 60 meV of the bosonic function (bottom-right inset of Fig. 5). The result (red 
solid line) is reported in Figure 5 and can be rationalized in terms of an impulsive decrease of 
the coupling and consequent narrowing of the Drude peak, that is, the opposite effect to the 
transient heating of bosons. 

In conclusion, we have shown that the use of non-equilibrium optical spectroscopy can 
be extended to doping and temperature regimes in which a pseudogap opens and electronic 
correlations strongly affect the electronic properties of the system. 

6. Conclusions and Perspectives 

Non-equilibrium optical spectroscopy is important and effective for investigating the physics of 
the cuprate superconductors and, more in general, of strongly correlated electron materials. 
Adding the temporal degree of freedom to the frequency resolution of equilibrium optical 
spectroscopy, allows disentangling the electronic and phononic contributions to the Eliashberg 
Bosonic Function I1(J7) in a prototypical cuprate, i.e. optimally-doped Y-Bi2212. Our results 
demonstrate that, in principle, the strength of the coupling of quasiparticles to bosons of 
electronic origin can account for the high-critical temperature of the system (T c =96 K). This 
result supports a description of high-temperature superconductivity in cuprates in terms of an 
attractive interaction, mediated by virtual bosonic excitations of electronic origin, where the 
most prominent candidates are spin fluctuations and loop currents. Furthermore, the analysis of 
the data obtained from non-equilibrium optical spectroscopy relies on a formalism (the Extended 
Drude model) that can be easily extended to more complex cases, e.g., systems with a gap in 
the density of states. This is very promising for further experiments aimed at investigating the 
elusive nature of the region of the cuprate phase diagram, in which a "pseudogap" in the density of 
states is present. Finally, the techniques described here can be easily extended to other complex 
materials, such as different classes of transition-metal oxides and iron-based superconductors. 

These new experimental and theoretical tools, developed to investigate the electronic 
properties of correlated materials out of equilibrium, open intriguing perspectives for the next- 
years materials science. The modeling of the interaction of ultrashort light pulses with correlated 
materials is still at its infancy. The use of the "effective" temperatures and of the extended Drude 
models to describe the non-equilibrium dynamics on the sub-ps timescale has been proved to be 
a very effective tool to carry out a quantitative analysis of the data and single out the different 
contributions to the pump-induced modification of the dielectric function. Nonetheless, these 
approximations cannot account for other important processes that play an important role in the 
femtosecond dynamics of strongly correlated materials. In particular: 

-the origin of the interplay between high-energy Mott-like excitations at 1.5-2 eV and the 
onset of high-temperature superconductivity, recently demonstrated by non-equilibrium optical 
spectroscopy [2|, is still an open question directly related to the nature of the superconducting 
phase in cuprates. 

-the k-space distribution of the nascent population, after interaction with a femtosecond light 
pulse, can have strongly non-thermal features, such as an excess of excitations at the antinodes 
[H] . This feature could provide the access to novel non-thermal phases [121 031 0H 051 06] , in 
which superconductivity can be controlled by means of light pulses. 

The joint effort of emerging theoretical tools, like non-equilibrium DMFT |47| 08] and 



time-dependent Gutzwiller approach [4£iJ ED] and time-resolved experiments probing different 
physical properties, like time-resolved optical and photoemission spectroscopies and time- 
resolved electron diffraction, will provide the key to understand the role of the electronic 
correlations in controlling the ultrafast electronic properties in unconventional superconductors 
and transition-metal oxides. 
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